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Abstract 

The Sculptor and Carina Dwarf spheroidal galaxies were observed with the 
H.E.S.S. Cherenkov telescope array between January 2008 and December 
2009. The data sets consist of a total of 11.8 and 14.8 hours of high quality 
data, respectively. No gamma-ray signal was detected at the nominal po- 
sitions of these galaxies above 220 GeV and 320 GeV, respectively. Upper 
limits on the gamma-ray fluxes at 95% C.L. assuming two forms for the spec- 
tral energy distribution (a power law shape and one derived from dark matter 
annihilation) are obtained at the level of 10^^^ to 10~^^ cm~^s~^ in the TeV 
range. Constraints on the velocity weighted dark matter particle annihila- 
tion cross section for both Sculptor and Carina dwarf galaxies range from 
(av) ~ 10~^^ cm^s~^ down to (av) ~ 10~^^ cm'^s"^ depending on the dark 
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matter halo model used. Possible enhancements of the gamma-ray flux are 
studied: the Sommerfeld effect, which is found to exclude some dark matter 
particle masses, the internal Bremsstrahlung and clumps in the dark-matter 
halo distributions. 

Keywords: Gamma-rays : observations - Dwarf Spheroidal galaxy. Dark 
Matter, Sculptor, Carina, Sommerfeld enhancement 
PACS: 98.70.Rz, 98.56.Wm, 95.35.+d 



1. Introduction 

The Dwarf Spheroidal Galaxies (dSphs) of the Local Group are the most 
common satellites of the Milky Way and assumed to be gravitationally bound 
dominantly by Dark Matter (DM). Although the predicted very high energy 
(VHE, E > 100 GeV) gamma ray flux from DM annihilation from dwarf 
galaxies is smaller compared to the expected DM annihilation gamma ray 
flux from denser regions of DM such as the Galactic Center, these galax- 
ies are promising targets for searches for gamma rays from DM annihilation 
since they are environments with a favorably low astrophysical gamma ray 
background. The galaxies themselves are expected to contain no astrophys- 
ical gamma ray sources since no recent star formation activity gives rise to 
VHE gamma rays (supernova remnants, pulsar wind nebula, etc.) and little 
or no gas acting as target material for cosmic rays has been measured 
Additionally their position at high galactic latitude is well separated from 
the numerous sources of VHE gamma rays in the Galactic plane. Also at 
such high altitudes no detectable contamination due to diffuse VHE gamma 
ray emission is expected, which in any case has shown up so far only from 
the Galactic Center region j2|. 

The H.E.S.S. {High Energy Stereoscopic System) array of Cherenkov tele- 
scopes has already observed dSphs and the collaboration has published re- 
sults on the Sagittarius dSph jH, 01 and the overdensity Canis Major jsj. 
In January 2008 H.E.S.S. launched observation campaigns on the Sculptor 
(^2000 = 01^02'^19"2, ^2000 = -33° 33' 00"0, distance ~ 79 kpc, [l|) and 
Carina (asooo = 06Ml"^36"0, ^2000 = -50° 58' 12"0, distance ~ 101 kpc, Q) 
dSphs, which are among the most luminous dSphs near the Milky Way. The 
Sculptor dSph was discovered in 1938 j6|, and was the flrst example of this 
type of galaxy in the vicinity of the Milky Way. The Carina dSph was dis- 
covered in 1977 01 . The best estimates of the orbits of the two dSphs show 
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that Carina is likely to be more tidally disrupted than Sculptor js], [oj , leading 
to higher uncertainties for the DM content of the Carina dSph than of the 
Sculptor. 

This paper presents the first results of a search for VHE gamma rays 
from DM annihilation from the Sculptor and Carina dSphs. It is organized 
as follows: in Section H] the analysis of the data is presented, from which up- 
per limits on the gamma ray flux are extracted assuming power-law spectra 
(Section [3]). In Section |4] flux limits for DM annihilation spectra are derived 
for both dSphs. Then the possibility of giving constraints on the DM particle 
properties is discussed (Section |5]), by assuming several DM galactic halo pro- 
files of the dSphs, and considering various possibilities for the DM candidate 
particle, in particular those that could give rise to an enhancement to the 
gamma ray annihilation flux. The results obtained are discussed in Section |6l 



The Fermi-LAT [Fermi-Large Area Telescope) collaboration 10| has recently 



published a search for gamma ray emission from the Sculptor dSph in the 



energy range from 100 MeV to 100 GeV, the results from Fermi-LAT [11 
will also be addressed in Section O 

2. H.E.S.S. observations and analysis 

2.1. The H.E.S.S. instrument 

The H.E.S.S. experiment is an array of four identical imaging atmospheric 
Cherenkov telescopes, observing VHE gamma rays. H.E.S.S. is located in 
the Khomas Highland of Namibia (23°16'18" South, 16°30'00" East) at an 
altitude of 1800 m above sea level. Each telescope has an optical reflector 
consisting of 382 round facets of 60 cm diameter each, yielding a total mirror 



area of 107 m^ 12|. The Cherenkov light, emitted by charged particles in 
the electromagnetic showers initiated by primary gamma rays, is focused 
on cameras equipped with 960 photomultiplier tubes, each one subtending 



a field-of-view of 0.16° [13|. The large field-of-view (~5°) permits survey 
coverage in a single pointing. The direction and the energy of the primary 
gamma rays are reconstructed by the stereoscopic technique [l3 ]. 

2.2. Observations and data analysis 

The observations of the Sculptor and Carina dSphs were conducted be- 
tween January 2008 and December 2009. They were performed in wobble 
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mode 15|, i.e. with the target typically offset by 0.7° to 1.1° from the point- 



ing direction, allowing simultaneous background estimation in the same field- 
of-view. The data used for the analysis were taken at average zenith angles 
of ~14° and ~34° for the Sculptor and Carina dSphs, respectively, leading 
to different effective energy thresholds. A minimal energy (-Emin) is defined 
as the energy for which the acceptance of the instrument reaches 20% of its 
maximum value, which for the Sculptor and Carina gives -Emin ~ 220 GeV 
and -Emin ~ 320 GeV, respectively. The data sets used for the analysis in- 



clude only the observations that meet standard quality control criteria [15 
based on the weather and data acquisition conditions. The total data set 
amounts to 11.8 h for Sculptor and 14.8 h for Carina of live time after the 
quality selection. These parameters are summarized in Table [T] 



dSph Galaxy 


Sculptor 


Carina 


Observation Period 


2008 Oct - 2008 Nov 


2008 Jan - 2009 Dec 


Live time (h) 


11.8 


14.8 


NoN 


117 


86 


Nqff 


2283 


1858 


a 


0.04 


0.05 


Significance 


I.Oct 


-1.4a 


Ar957o U.L. 
^^7, tot 


32.4 


8.6 


^min (GeV) 


220 


320 




19.2 


10.2 



Table 1: H.E.S.S. observation characteristics and upper limits on the ob- 
served number of gamma rays for the Sculptor and Carina dSphs. Nqn and 
Nqff are the number of gamma-ray candidate events in the signal region 
and in the background region, respectively, a is defined as the ratio of the 
on-source area to the off-source area. The significance of the excess in the 
signal region is calculated for the given 

NoN, NoFF and a. N^%^-^- is the 
95% confidence level upper limits on the total observed numbers of gamma- 
rays, and N^^'^'''^-^-{E.y > i?min) is the 95% confidence level upper limits on 
the observed numbers of gamma-rays above the given minimal energy -Bmin 
for each dSph. The minimal energy -Emin is defined as the energy for which 
the acceptance of the instrument reaches 20% of its maximum value. 



The data are analyzed using the model analysis from [16|], with standard 
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cuts. This method is based on an accurate pixel per pixel comparison of the 
observed intensity with a pre-calculated semi-analytical model of showers. 
It provides an energy resolution of ~10% with an angular resolution at 2a, 
defined as the 95% containment radius, of the order of magnitude of 0.1° 
in most of the energy range. The background was determined by the ring- 



background technique [17[, calculating the background for each position in 
the field-of-view using the background rate contained in a ring around the 
target. 

No significant gamma ray excess was found above the estimated back- 
grounds at the nominal positions of Sculptor and Carina dSphs, as seen in 
Figure m The significance of the excess in the 0.1° radius integration area 
for Sculptor and Carina are -1-1.0 a and —1.4 a, respectiveljU. This allows 
one to set 95% confidence level upper limits on the total observed numbers 
of gamma-rays N^^'^^^'^' , using the number of gamma-ray candidate events 
in the signal region Nqn, in the background region Nqff, and the ratio a 
of the on-source area to the off-source area. The limits have also been com- 
puted above the given minimal energy E^in for each dSph. The iV^^^^.L. 
and iVf ^"C.L^^^ > ^^.j obtained using the method described in Feldman 



& Cousins 19| are shown in Tabled! 



3. Flux upper limits with power-law gamma ray spectra 

To obtain flux upper limits for standard astrophysical sources, power-law 
photon flux spectra of index F are assumed. 

The resulting 95% confidence level upper limits on the gamma ray fluxes 
^95% c.L. j^-g^ > Ejnin) are derived through the knowledge of the acceptance of 
the instrument for each energy bin, given the parameters of the observation 
data sets, e.g. the mean zenith angle, the optical efficiency, the mean offset 



angle and the cuts used in the analysis [15[. The index F was varied between 



■^A sensitivity of 20.1 events is reported for Carina, following the suggestion in Feldman 
& Cousins [l^ in the case where fewer events than the expected background are measured. 
The sensitivity is defined as the average upper limit (at 95% C.L.) obtained with the 
expected background and no true excess signal. 
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1.8 and 2.4, which correspond to standard spectral indexes for astrophysical 



sources 



20| . The results are summarized in the Table [2J 



Galaxy 


$95%c.L.(E^ > E^in) (10-13 cm-^s-i) 


Sculptor (Einin = 220 GeV) 
Carina (E^in = 320 GeV) 


5.1 - 6.2 
1.6 - 2.0 



Table 2: The 95% confidence level upper limits on the gamma ray fiuxes 
above the minimal energy -Emin, given in units of 10-^^ cm-^s'^, for a power- 
law model with indices between F = 1.8 and F = 2.4. The lower values of 
the upper limits are found for the index F = 1.8. 



4. Flux upper limits with Dark Matter annihilation spectra 

4.I. Theoretical framework 

Theories beyond the Standard Model (SM) of particle physics propose 
several candidates to account for the DM inferred from the gravitational 
interaction. Self-annihilation of DM particles into ordinary matter will in 
general produce gamma rays as sub-products, which may be detectable with 
current instruments. 



For instance, some supersymmetric extensions of the SM [21[ predict the 



neutralino as the lightest stable supersymetric particle, which is a good candi- 



date for DM [2l|. Several parametrizations for the continuum contribution to 



the total neutralino self-annihilation gamma-spectrum are used here [22|, |23 
The latter arises from the hadronisation and decay of the cascading prod- 
ucts. Universal extra dimensional extensions of the SM provide also suitable 
DM candidates. In the Kaluza-Klein (KK) models the first KK mode of the 



hypercharge gauge boson B^^^ is a DM particle candidate 24 . 

The gamma ray fiux from the annihilations of DM particles of mass jtidm 
in a DM halo is given by a particle physics term times an astrophysics term: 

d$^(Afi,^^) 1 {av) dN. 



dE^ 8tt m'^j^ 

Particle Physics 

where the astrophysical factor is defined as 



X J{An)An, (2) 

Astrophysics 



J{An) = ^ [ dQ [ py{s)]ds. (3) 

t/ j^S^ t-* los 
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In equation [3] the squared density of DM (p^) is integrated along the hne 
of sight (los) and over the sohd angle AQ. The solid angle is fixed as the 
angular resolution of the telescope for a point-like source search. For the 
H.E.S.S. experiment AO = 10~^ sr. The particle physics term contains the 
DM particle mass, mDM, the velocity- weighted annihilation cross section, 
{(jv), and the differential gamma ray spectrum from all final states weighted 
by their corresponding branching ratios, dN^/dE^. 

The number of detected gamma rays above a minimal energy -Emin is 
given by 

N^{E > E^^) = Tobs / A,s{E^) dE^ , (4) 

where Tobs is the observation time, and Ac{[{E^) is the effective area of the 
detector as a function of the gamma ray energy, the zenith angle, the offset of 
the source from the pointing direction and the selection cuts. The differential 
flux d^j/dEj, in the case of the DM annihilation, is given by equation |2l 

4-2. Flux upper limits with neutralino annihilation spectra 

Before modelling the shape of the DM halos and addressing the astro- 
physics factor J{AQ)AQ, upper limits on the flux at 95% CL. for different 
annihilation spectra are considered as a function of the DM particle mass. 
The number of gamma rays above a minimal energy can be related to the 
integral flux by 



[E > E^^) ^ X $^ [E> E^^) . 



(5) 

The annihilation spectrum depends on the nature of the final states. A 
parametrization using the the average of the WW and ZZ final states was 



taken from Bergstrom et al. 22|, which will be used here and in Section 



For annihilation into bb and r+r~ PYTHIA 6.225 [25[ was used to compute 
the spectra. The upper limit on the flux at 95% C.L. is related to the upper 
limit on the number of observed gamma rays at 95% C.L. by replacing the 
jy95%c.l. f^Q^ Section O (Table OP) in equation O 

Figure |2] shows the calculated upper limit on the flux for both Sculptor 
and Carina dSphs, which depends on the assumed spectrum and hence on 
the mass of the neutralino. The results obtained by Fermi-LAT for the 
Sculptor dSph and energies > 100 MeV are also plotted. As can be seen. 
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for high neutrahno masses (> 500 GeV) H.E.S.S. is more sensitive. The flux 
sensitivity is quahtatively driven by the product of the acceptance Aeff{E^) 
times the observation time Tobs- Using the acceptances of about Aes ~ 10^ m^ 
for H.E.S.S. and of a few m^ for Fermi- L AT, and observation times of about 
~12 hours for H.E.S.S. and ~11 months for Fermi-LAT, the ratio between 
their sensitivities for a given DM mass yields a better sensitivity for H.E.S.S. 
by a factor of a few hundred, for masses well above the H.E.S.S. threshold. 

5. Exclusion limits on the Dark Matter annihilation cross section 

5.1. Astrophysical factor calculation 

In order to calculate the exclusion limits on the DM annihilation cross 
section, one needs to model the density distribution of DM in the observed 
target that will be used in the astrophysical factor J calculation. Two hy- 
potheses for spherical DM halo profiles are used for Sculptor and Carina: a 



pseudo-isothermal profile 26|, and the Navarro, Frenk, and White (NEW) 



profile [27 



The pseudo-isothermal sphere case belongs to the cored profile class. It 
is characterized by the mass density profile 

Piso(r) = po ^^,^^,^ , (6) 

where is the core radius and po is a normalization factor. This profile is 
completely defined by po and Tc- 

The NEW profile is a cusped profile, with the mass density given by 

Pnfw('^ ) = -ri~Tn~7 — ' (''') 

(r/rs)(l + r/r^y 

where is a scale radius, p° the present critical density of the Universe and 
6c a characteristic over-density, defined by 6c = c^g{c)/3, with = 200. 
Here c = R^/r^ is the concentration parameter of the halo, where is the 
virial radius, and g{c) = (ln(l + c) — c/{l + c))~^. In this case the profile 
is completely defined by the concentration parameter c and the virial mass 
My, expressed as 

= y A.p°i?^ . (8) 



= 2.775 X 10" h"^ MoMpc"^ = 1.053 x 10"^ /^2 (QeV/c^) cm^^ 
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Given a hypothesis about the gravitational potential of a galaxy, i.e. its DM 
halo profile, and about the velocity dispersion anisotropy of its stars, one can 
obtain the theoretical line-of-sight velocity dispersion from the Jeans equa- 
tion. The methodology used to find each set of parameters for each halo 
profile used in the literature cited below consists of comparing the observed 
line-of-sight velocity dispersion of stars, o"ios, for each distance bin with that 
predicted for the various models. Minimizing the between the model and 
the data provides the best set of parameters. o"ios is also dependent on the 
stars velocity dispersion anisotropy /3{r) (see e.g. equation 7 of 29[). The 
latter is taken into account in the fitting procedure and different parametriza- 
tions were used. 



Sculptor: The DM halo mass content of Sculptor dSph as well as its 
profile parameters are estimated in a recent paper 



30 and thesis 31 



us- 



ing two RGB [Red Giant Branch) stars populations to partially break the 
Jeans degeneracy in the DM halo modelling 32|. The last gravitationally 
bound star was measured at rjast = 1-8 kpc, which gives an estimate of 
the tidal radius j3l[ and is used in the integration to obtain J(Af2). Two 
hypotheses for the /3(r) profile were explored in [3o[: a radially constant ve- 
locity dispersion anisotropy, and a Osipkov-Merritt (OM) velocity dispersion 



anisotropy (33|, |3J]. Using the parameters extracted from [3l|, the astro- 
physical factor J is computed for eight different Dark Matter halos. The 
parameters as well as the astrophysical factor are summarized in the Ta- 
ble [3l The large variety of investigated DM halos allows one to encompass 
the astrophysical uncertainties induced by the modelling. 



Carina: Here, the best fit parameters of each DM halo profile were 
given for a NFW profile in 35| and for the pseudo-isothermal profile in 36 



The parameters for the NFW profile obtained from [35| are the virial mass 
My = 2 X 10^ Mq and the star's velocity dispersion anisotropy /3(r) = —0.5. 
Replacing in equation [HI the virial radius is found to be = 12 kpc. 
Using the relationship between and the concentration parameter c found 
in [33] (see also j38|) yields c ~ 22. Finally rg = 5.35 x 10~^ kpc is obtained 
using rg = c/R^. The parameters for a pseudo-isothermal profile with an 
isotropic velocity dispersion (/3(r) = 0) and consistent with the kinematic 



39| . see also Figure 4 of (36 



data were obtained from 

set arbitrarily to rt = 2.0 kpc, which is a conservative value 



The tidal radius is 
4l|, and it 



40 



is used in the integration to obtain J[/SD). The parameters of the DM halo 
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Sculptor dSph 


Iso Profile 


rc 
(kpc) 


' last 

(108 Mq) 


Po 

(10^ MQkpc-3) 


J 

(1023 GeV^cm-S) 


P = const 


0.05 


1.2 


221 


2.98 


0.5 


3.3 


9.13 


0.27 


P = Pom 


0.05 


1.3 


240 


3.49 


0.5 


3.4 


9.40 


0.29 


NFW Profile 


c 


(109 Mo) 


(kpc) 


J 

(1023 GeV2cm-5) 


P = const 


20 


1.9 


1.26 


2.75 


35 


0.59 


0.48 


5.20 


P = Pom 


20 


2.2 


1.32 


3.29 


35 


0.68 


0.51 


6.24 



Table 3: Structural parameters of the eight best fits (30| in the case of a 
pseudo-isothermal and NFW DM halo profiles for the Sculptor dSph, as well 
as the corresponding values of the astrophysical factor J, for two hypotheses 
of the velocity anisotropy profile P{r), two core radii (pseudo-isothermal) and 
two concentration parameters (NFW). 



profiles as well as the astrophysical factor J are summarized in Table HI 



Carina dSph 


Iso Profile 

p = o 


rc 
(kpc) 


rt 
(kpc) 


Po 

(10^ M0kpc-3) 


J 

(1022 GeV2cm-5) 


0.22 


2.0 


1.36 


2.01 


NFW Profile 
P = -0.5 


c 


(109 Mo) 


T 

(kpc) 


J 

(1022 GeV2cm-5) 


22 


0.20 


0.54 


4.37 



Table 4: Structural parameters of the two best fits 35|, l36|] in the case of a 
pseudo-isothermal and NFW DM halo profile for the Carina dSph, as well as 
the corresponding value of the astrophysical factor J (see equation [3]) . 
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5.2. Generic case for exclusion limits 

The 95% CL. upper limit on the velocity-weighted annihilation cross 
section as function of the DM particle mass for a given halo profile is given 
by 

where the parametrization of dNy/dE^ is taken from j22[ for typical neu- 
tralino self-annihilation into WW and ZZ pairs, and calculated from 2J| for 
Kaluza-Klein B^^^ self-annihilation. The exclusion curves for the neutralino 
case are plotted for the Sculptor and Carina dSphs in Figure [3] referring to the 
halo profiles given in the Table [3] for Sculptor, and in the Table H] for Carina, 
respectively. The Fermi-LAT exclusion limit for Sculptor is added extend- 



ing up to 1 TeV pjj, which is based on a NFW profile with = 0.9 kpc 



and ps = X Pc = 3.7 X 10 Mokpc , and a neutralino parametrization 
with only bb in the final state. Using the parameters and ps from Fermi- 



LAT paper [ll|, the astrophysical factor with the H.E.S.S. solid angle is 
J = 1.33 X 10^^ GeV^cm"^. The resulting exclusion limits are plotted (pink 
dashed line). 

Below ~ 1 TeV, the Fermi-LAT results provide stronger limits than the 
H.E.S.S. results. In comparison with the flux sensitivity (Section 14. 2p . the 
{av) upper limits also take into account the predicted integrated number of 
gamma rays in the instrument energy range. The predicted number of gamma 
rays per annihilation event in the Fermi-LAT energy range is about 10^ times 
higher than the one in the H.E.S.S. energy range. This implies a Fermi-LAT 
limit which is of the order of 10 times better than the one for H.E.S.S., 
despite the latter's stronger flux sensitivity. Fermi-LAT and H.E.S.S. give 
complementary limits on (av) in the 10 GeV - 100 TeV mass range. Since 
the given limits on the cross section are a few orders of magnitude above the 
theoretical expectations for DM particles in the Minimal Supersymmetric 
extension of the Standard Model (MSSM), those expectations are not shown 
in the figures. 

The Figure m shows the exclusion limits of {av) in the case of the Kaluza- 
Klein DM particle B^^\ The limits are plotted for the Sculptor dSph referring 
to the halo profiles given in the TableO In the TeV range the 95% C.L. upper 
limit on {av) reaches 10~^^ cm^s"^. 
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5. 3. Enhancement effects for the exclusion limits 

Three cases that can modify the exclusion hmits are considered: two par- 
ticle physics effects, namely the Sommerfeld enhancement and the Internal 
Bremsstrahlung (IB) from the DM annihilation, and an astrophysical effect 
due to the mass distribution of dark-matter sub-halos. 

5.3.1. The Sommerfeld enhancement 

Here two new assumptions were made for the Sculptor's and Carina's 
DM halo composition. The first assumption is to assume the DM particle 
annihilates to a W boson, which is the case when the neutralino is a pure 
wino 42|. The second assumption is to assume that the DM mean velocity 
inside the halo is the same as for the stars (the mean velocity dispersion of 
the stars is ay ~ 10.0 km/s for Sculptor and cry ~ 7.5 km/s for Carina). The 
later assumption is plausible due to the relaxation time scales in the dwarf 
galaxies. 

In this class of objects, the relative velocity between the DM particles 
may be sufficiently low so that the Sommerfeld effect can substantially boost 



the annihilation cross section |43|, since it is particulary effective in the very 
low-velocity regime. This non-relativistic effect arises when two DM particles 
interact in an attractive potential. The actual velocity-weighted annihilation 
cross section of the neutralino can then be enhanced by a factor S defined as 

{av) = S {av), , (10) 

where the value of S depends on the DM particle mass and relative velocity. 
A wino would interact with the attractive potential created by the Z gauge 
boson through the weak force before annihilation occurs, which would give 
rise to an enhancement. The value of this enhancement was numerically 



calculated as done in j43| and then used to improve the 95% C.L. upper limit 
on the velocity-weighted annihilation cross section, (av) /S* as a function of 
the DM particle mass. The effect of this enhancement is shown in Figure [5] 
for Sculptor (top) and Carina (bottom), and for two particular cases of the 



halo profile models. The predicted {av)^ for a pure wino j42| as well as 



the typical annihilation cross section for a thermally produced DM {{(tv)q 



10~^^ cm^s~^ 2l|) are also plotted. Some specific wino masses can be 
excluded at the level of {a'v)^ ~ 10~^^ cm^s~^. 

A more general case of the enhancement by the Sommerfeld effect can 
be treated by changing the exchanged boson mass and the coupling con- 
stant accordingly. The position of the resonances is qualitatively driven by 
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equation 5 of [43|. Increasing the boson mass shifts the resonance to higher 
DM masses since the the weak couphng constant is only a slowly varying 
function of the boson mass. Simultaneously, once the relative velocity and 
the DM particle mass are fixed, the value of the enhancement close to the 



resonance grows roughly linearly with the boson mass |43| . Examples of the 



enhancement for 0.1 to 1.0 GeV exchange boson masses are shown in (44 



5.3.2. Internal Bremsstrahlung 

In some specific regions of the MSSM parameter space, the electromag- 
netic radiative correction to the main annihilation channels into charged 
particles can give a significant enhancement to the expected gamma ray flux 
in the observed environment due to internal Bremsstrahlung (IB) |45l . |46 |. 
In the stau co- annihilation region of the minimal supergravity (mSUGRA) 
models, for instance, the wino annihilation spectrum would receive a consid- 
erable contribution from the internal Bremsstrahlung. 

This contribution to the annihilation spectrum was computed using the 



parametrization of [46[ for all the wino masses in the H.E.S.S. energy range. 
The enhancement effect on the 95% C.L. upper limit on the velocity-weighted 
annihilation cross section is shown in Figure El The joint enhancement due 
to the Sommerfeld effect and IB is also plotted. The effect of the IB is only 
significant in the exclusion limits for the low mass DM particle regime. 

5.3.3. Enhancement from dark-matter sub-halos 

Astrophysical effects may also modify the exclusion limits. Numerical 
simulations of galactic halos predict a population of subhalos that could con- 
tribute to the overall astrophysical factor in equation [31 Using the procedures 



given in |47[, the contribution to the astrophysical factor by the DM sub-halos 
population is estimated. An enhancement of the astrophysical factor is found 
to be of a few percent, which is too small to significantly affect the exclusion 
limits presented. 



6. Summary 

Both Sculptor and Carina dSphs are well-studied in multiple wavelengths, 
providing reasonable measurements of the profile of the DM in their halos. 
Recent VHE gamma ray observations from H.E.S.S. of both of these objects 
provide new insight into the dark matter within them. While no positive 
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dark-matter annihilation signal has been detected, the observations provide 
constraining limits on dark matter parameters. 

Constraints have been obtained for the velocity weighted annihilation 
cross section {av) as a function of the mass for neutralino and KK DM 
particles. Concerning Sculptor dSph upper limits on {av) have been ob- 
tained in the range of ~10~^^ cm^s"^ to ~10~^^ cm^s~^ for neutralinos and 
~10~^^ cm'^s"^ to ~10~^'^ cm^s~^ for KK particles. From the different profile 
parameters investigated, much better limits are obtained for a NFW profile 
with a strong concentration parameter c = 35 when compared to the limits 
obtained for a isothermal profile with a large core radius of Tc = 0.5 kpc. Also 
for the Carina Sph better limits for (av) of neutralinos have been obtained 
for a NFW profile. 

The DM halo model induces systematic uncertainties in the exclusion 
limits: the value of the astrophysical factor can vary over one order of mag- 
nitude for a given halo profile in the case of Sculptor. The results presented 
show that the DM particle models that could satisfy WMAP constraints on 



the Cold Dark Matter relic density [48| cannot be tested. H.E.S.S. limits 



are comparable to the limits reported by MAGIC |49| and VERITAS [50 
but weaker than those obtained by Fermi- L AT ll| in the GeV mass range. 
Nevertheless, they are complementary to the Fermi-LAT limits in the TeV 
range. Searches through neutrinos provide competitive limits. In the WW 



channel IceCube constraints on (erf) lie at the level of ~10~^^ cm^s~^ 51 
Finally assuming a resonance effect in the Sommerfeld enhancement, some 
specific wino masses can be excluded, and the first experimental constraints 
have been obtained on the Sommerfeld effect using H.E.S.S. data and DM 
annihilation spectra. 
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Figure 1: Oversampled significance maps (with an integration region of 0.1°, 
wfiite circle) in equatorial coordinates calculated according to the Li & Ma 
method [isj in the directions of Sculptor (top left) and Carina (bottom left) 
dSphs. Distributions of the significance in the maps for the Sculptor (top 
right) and the Carina (bottom right) dSphs. The solid line is a Gaussian 
fitted to the data. The significance distribution is well described by a normal 
Gaussian. No significant excess is seen at the target position. 
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Figure 2: (top) Upper limits on fluxes at 95% CL. for the Carina (dashed 
lines) and Sculptor (solid lines) dSphs obtained with H.E.S.S., for E > 
320 GeV and E > 220 GeV, respectively, as function of the DM mass. The 
parametrization of the annihilation spectrum using the the average of the 
WW and ZZ final states is extracted from Bergstrom et al. 2^. Fermi-LAT 
results for Sculptor with E > 100 MeV are also plotted for t~^t~ and bb 



final states [ll|. (bottom) A zoomed view on the upper limit from H.E.S.S. 



using in addition various parametrizations for the annihilation spectrum [25 
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Figure 3: Upper limit at 95% CL. of (av) as function of the DM particle mass 
for different DM halos for Sculptor {top) and Carina [bottom). For the NFW 
halo profile of Sculptor two concentration parameters are used: 20 and 35. 
For the pseudo-isothermal halo profile two core radii are used: 0.05 kpc and 
0.5 kpc. Two hypotheses on the velocity anisotropy parameter are also stud- 
ied: a constant (solid lines) and an Osipkov-Merritt (dashed lines) anisotropy. 
The velocity anisotropy and the concentration parameters are given in brack- 
ets for the NFW profile. The velocity anisotropy and the core radius are given 
in brackets for the pseudo-isothermal profile. The Fermi-LAT limits ll| for 
a NFW profile are also plotted as well as the H.E.S.S. limits for this NFW 
profile (rs = 0.9 kpc and = 3.7 x 10^ Mgkpc'^). For Carina both the NFW 
halo profile and the pseudo-isothermal halo profile are plotted (see text for 
parameters) . 
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Figure 4: Upper limit at 95% CL. of {av) as function of the Kaluza-Klein 
DM particle mass for all DM profiles of Sculptor in Table [31 The veloc- 
ity anisotropy and the concentration parameters are given in brackets for 
the NFW profile. The velocity anisotropy and the core radius are given in 
brackets for the pseudo- isothermal profile. 



23 





Figure 5: Upper limit at 95% CL. on [av) / S as function of the DM parti- 
cle mass enhanced by the Sommerfeld effect (see text for more details) for 
Sculptor {top) and Carina (bottom). The NFW halo profile as well as the 
pseudo-isothermal profile are used. The predicted {(tv)q for a pure wino jiij 
(solid black line) as well as the typical cross section for a thermally produced 
DM (dashed red area) are also plotted. 
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Figure 6: Upper limit at 95% CL. on {av) / S as function of the DM particle 
mass enhanced by the Sommerfeld effect and the internal Bremsstrahlung 
(see text for more details) for a NFW profile of Sculptor. The predicted 
((Tw)q for a pure wino (solid black line) as well as the typical cross section for 
a thermally produced DM (dashed red area) are also plotted. 
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